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Abstract 

For a general Multidimensional Levy process (satisfying some moment conditions), we introduce the Multidimen- 
sional power jump processes and the related Multidimensional Teugels martingales. Furthermore, we orthogonalize 
the Multidimensional Teugels martingales by applying Gram-Schmidt process. We give a chaotic representation for 
every square integral random variable in terms of these orthogonalized Multidimensional Teugels martingales. The 
predictable representation with respect to the same set of Multidimensional orthogonalized martingales of square in- 
tegrable random variables and of square integrable martingales is an easy consequence of the chaotic representation. 
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1. Introduction 

The chaotic representation property(CRP) has been studied by Emery(1989) for normal martingales, that is, for 
martingales X such that < X, X >,- ct, for some constant c > 0. This property says that any square integrable ran- 
dom variable measurable with respect to X can be expressed as an orthogonal sum of multiple stochastic integrals 
with respect to X. It is known (see for example Dellacherie et al., 1992, p. 207 and Dermoune, 1990), that the only 
normal martingales X, with the CRP, or even the weaker predictable representation property (PRP), which are also 
Levy processes are the Brownian motion and the compensated Poisson process. David Nualart and Wim Schoutens 
(2000) study the chaotic representation property for one-dimensional Levy process, in terms of a suitable orthogonal 
sequence of martingales where these martingales are obtained as the orthogonalization of the compensated power 
jump processes of the Levy process. Furthermore, Nualart and Schoutens (2001) used their martingale representation 
result to establish the existence and uniqueness of solutions for BSDE's driven by a Levy process of the kind consid- 
ered in Nualart and Schoutens (2000). In addition, Corcuera, Nualart and Schoutens (2005) applied this martingale 
representation result to study the completion of a Levy market by power-jump assets. 

In the past twenty years, there is already a growing interest for multidimensional Levy Processes. Some con- 
cepts and basic properties about multidimensional Levy Processes were summarized in Sato (1999). Applications 
of multidimensional Levy Processes to analyzing biomolecular (DNA and protein) data and one-server light traffic 
queues were explored by Dembo, Karlin and Zeittouni (1994). A small deviations property of multidimensional Levy 
Processes were discussed by Simon (2003). In finance research, practically all financial applications require a mul- 
tivariate model with dependence between components: examples are basket option pricing, portfolio optimization, 
simulation of risk scenarios for portfolios. In most of these applications, jumps in the price process must be taken into 
account. Cont and Tankov (2004) systematically investigated these problems in multidimensional Levy market. In ad- 
dition, the optimal portfolios in multidimensional Levy market is discussed by Emmer and Kliippelberg (2004). Some 
simulation approaches for multivariate Levy processes are also investigated in Cohen and Rosihski (2007). Levy 
copulas was also suggested by Kallsen and Tankov (2006) in order to characterize the dependence among components 
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of multidimensional Levy Processes. The SDEs driven by infinite-dimensional Levy processes was investigated by 
Meyer-Brandis and Proske (2010). 

The chaotic representation property is important for the research of Levy process, and multidimensional Levy pro- 
cesses obtain some applications in bioscience and finance, so it is significant to extend the result in univariate set-up 
obtained by Nualart and Schouten(2000) to the cases of multidimensional Levy processes. In this paper, following 
the research line of the paper in Nualart and Schouten(2000), we study the chaotic representation property for Multi- 
dimensional Levy processes, in terms of a suitable orthogonal sequence of Multidimensional martingales, assuming 
that the Levy measure has a finite Laplace transform outside the origin. These Multidimensional martingales are ob- 
tained as the orthogonalization of the Multidimensional compensated power jump processes of our Multidimensional 
Levy process. In Section 2, we introduce these Multidimensional compensated power jump processes and we trans- 
form them into a multivariate orthogonal sequence. Section 3 is devoted to prove the chaos representation property 
from which a predictable representation is deduced. Finally, in Section 4, we discuss some particular examples. 

2. Preliminary 

A Revalued stochastic process X = \X(f) = (Xi(t), X 2 (t), • • • , X n (t))', t > 0) defined in complete probability space 
(Q., & ', V) is called Levy process if X has stationary and independent increments and X(0) = 0. A Levy process 
possesses a cadlag modification (Protter, 1 990, Theorem 30,p.21) and we will always assume that we are using this 
cadlag version. If we let & t = % V JV, where % = cr{X(s), < s < t] is the natural filtration of X, and ,jV are 
the P-null sets of then t > 0} is a right continuous family of cr-fields (Protter, 1990,Theorem 31,p.22). We 
assume that & is generated by X. For an up-to-date and comprehensive account of Levy processes we refer the reader 
to Bertoin (1996) and Sato (1999). 

Let X be a Levy process and denote by 

X(t-)= lim X(s), f>0, 
s-»f,s<; 

the left limit process and by AX(t) = X(t) - X(t-) the jump size at time t . It is known that the law of X(t) is infinitely 
divisible with characteristic function of the form 

E [exp(W ■ X(t))] = (0(0))' , = (6»i, 6» 2 , • • • , n ) e W 

where 0(0) is the characteristic function of X(l). The function ifi(0) = log(f>(0) is called the characteristic exponent 
and it satisfies the following famous Levy-Khintchine formula (Bertoin, 1996): 

0(0) = -1.0 ■ 26* + ia ■ + f (exp(i0 x)-l-i0- xl M <i) v(dx). 

where a, x e W, 2 is a symmetric nonnegative-definite n x n matrix, and v is a measure on M"\{o} with J (\\x\\ 2 A 
l)v(dx) < oo. The measure v is called the Levy measure of X. 

Throughout this paper, we will use the standard multi-index notation. We denote by No the set of nonnegative 
integers. A multi-index is usually denoted by p, p = (p\,pi, • ■ • ,p n ) £ Whenever p appears with subscript or 
superscript, it means a multi-index. In this spirit, for example, for a: = (x\ , ■ • ■ , x n ), a monomial in variables x\ , ■ ■ ■ ,x n 
is denoted by x p = x p ^ ■ ■ ■ x„" . In addition, we also define p\ = p\ ! • • ■ p n \ and \p\ = p\ + ••• + p„; and if p, q e Ng, 
then we define S p ^ q = 6 Puqi ■ ■ ■ 5 Pn ,q n - 

Hypothesis 1. We will suppose in the remaining of the paper that the Levy measure satisfies for some s > 0, and 
A > 0, 

I exp(A\\x\\)v(dx) < oo. 

J\X\>£ 

This implies that 

j x p v(dx) < oo. \p\ > 2 (1) 
2 



and that the characteristic function E [exp(i6 ■ X(t))] is analytic in a neighborhood of origin o. As a consequence, X(t) 
has moments of all orders and the polynomials are dense in L 2 (R", P o X(f) -1 ) for all t > 0. 

Professor Nualart proposed author to use the following transformations of X which will play an important role in 
our analysis. We introduce power jump monomial processes of the form 

x(t) (Pu-, P „) « £ (aXi(s)) Pi ■ ■ ■ (AX n (s)r% 

Q<s<t 

The number \p\ is called the total degree of X(t) p . Furthermore define 

Y (f)(pi.-Jx.) d £f xw (p,, -,p n ) _ E[X(f) (/ "' , "' p " ) ] = Xit)^'-'"'^ - m p t, 

the compensated power jump process of multi-index p = (p\,p2, • • • ,p n )- Under hypothesis 1, Y{t) (j ' l '' ,Pn) is a normal 
martingale, since for an integrable Levy process Z, the process \Z, - E[Z t ], t > 0} is a martingale. We call Y(t) (j?l, '" ,Pn ^ 
the Teugels martingale monomial of multi-index {p\, ■ ■ ■ ,p n ). 

Remark 1. In the case of a Poisson process, all power jump processes will be the same, and equal to the original 
Poisson process. In the case of a Brownian motion, all power jump processes of order strictly greater than one will be 
equal to zero. 

In the following we will introduce some concepts and basic properties of martingale polynomial. These new 
concepts and properties are totaly similar to those of polynomial in n real variables (cf. Dunkl and Xu (2001)). 

A martingale polynomial P in n Levy variables X = (X\,X2, ■ ■ ■ ,X n ) is a linear combination of martingale mono- 
mials, 

P(X) = Yj c p Y "> 
IpI>i 

where the coefficients c p are in the real numbers R. The degree of a martingale polynomial is defined as the highest 
total degree of its martingale monomials. We shall use the abbreviation II" to denote the collection of all martingale 
polynomials in X. We also denote the space of martingale polynomials of degree at most d by ITJ. A martingale 
polynomial is called homogeneous if all the monomials appearing in it have the same total degree. Denote the space 
of homogeneous polynomials of degree d e N in n variables by !PJ; that is 

K = \P:P(X)=Y j c p yA. 

{ \p\=d J 

Every polynomial in II" can be written as a linear combination of homogeneous martingale polynomials; for P e ITj, 

p{x) = f j Y J c P y. 

k=\ \p\=k 

Denote by r" d the dimension of P" d and it is well known that 

r» = dim<r d = [ d + n d ~ 1 ) and dimll" = ( d + n ) - 1. 

We denote by jV 1 the space of one dimensional square integrable martingales M such that sup E(M(f) 2 ) < <x>, and 

t 

M(0) = a.s. Notice that if M e ,yV 2 , then limE(M(f) 2 ) = E(M(oo) 2 ) < oo, and M(t) = E[M(oo)|^]. Thus, each 

t—>co 

M e .yV 1 can be identified with its terminal value M(oo). As in Protter(2005, p. 181), we say that two martingales 
M, N e are strongly orthogonal and we denote this by M x N, if and only if the product MN is a uniformly 
integrable martingale. As noted in Protter (2005,pl81), one can prove that M x N if only if [M,N] is a uniformly 
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integrable martingale. We say that two random vectors X,Y e L 2 (Q, are weakly orthogonal, X ± Y,if E[XY] = 0. 
Clearly, strong orthogonality implies weak orthogonality. 

In the theory about orthogonal polynomials of several variables, we can apply the Gram-Schmidt process to 
the monomials with respect to the usual inner product to produce a sequence of orthogonal polynomials of several 
variables. Some details about the technique and theory of orthogonal polynomials of several variables refer to Dunkl 
and Xu (2001). In this paper, we shall apply the standard Gram-Schmidt process with the graded lexicographical order 
to generate a biorthogonal basis {H p , p e N"), such that each H p (\p\ = d) is a linear combination of the Y q e Yl n d , with 
\q\ < \p\ and the leading coefficient equal to 1. We set 

W-i 

HP = Y? + J] c q Y" + Y u ^ c " Y ^ 

q<p,\q\=\p\ k=\ \q\=k 

where p = {p\, ■ ■ ■ ,p„], q = \qi, - ■ ■ ,q„] and < represent the relation of graded lexicographical order between two 
multi-indexes. 
We have that 

[HP, Y"](t) = 2 (AX^s))^ ■ ■ • (AX„( f ))"" + «" 

0<s<f 

+ 2 2 c p+q {AX,{s))p^ ■ ■ ■ (AX n (s)y^- + ££ c ei+q o- u ti {q=ei] . 

0<.s<(l<|p|<|p| i=l 7=1 

n 

Let m p+q = J n *f qi v{dx), then we have that 

i=l 

/ n n \ 

E[H p ,Y q ](t) = t m p+q + ^ Cp +q mp +q + ^^c e . +q o-ijtI {q=e]] 
V 1<\P\<\P\ (=1 7=1 

where c, denotes the n-dimensional unit vector with ;'th component equal to one. In conclusion, we have that, [H p , Y q ] 
is a martingale if only if we have that E[H P , Y q ](l) = 0. 

Consider two spaces: The first space S i is defined as follows 

Si = l^^ctipu--- ,p n )x p > ---4" +c Q + ^ c_i0'i,--- ,i n )x'l •••4"; 

[k=l \p\=k (ii,-,fc,)e{0,-l)»,|i|>-(B-l) 

d e {l,2,---},cj(pu--- ,p n ) e R,j = -1,0, - •• ,d;x t * 0, i= 1,2, •••,«;/ = (i u ■■ • ,/„)} 
which is endowed with the scalar product < •, • >i given by 

X+co f*+oo n 

• I P(x)Q(x)Y\x*v(dx) 
CO xJ — OO ■_ j 

n n 

+ Z Z cic2 ^W)=c 1 ^- 1 , ew=C2 ^-i ) - i = a, i, ■ • ■ , i). 

'=1 7=1 

Note that 



+ Z Z (T U I {p=e„q=ejh \P\ > L \q\ > L 



1=1 7=1 



Thus we can construct the other space 5 2 which is the space of all linear transformations of the Teugels martingale 
monomials of the multivariate Levy process, i.e. 

S 2 ={ 2 c d (p u --- ,p„W(t)(P--^+ 2 Cii-\(p\, - ■ ■ , p„)Y(tY p, '"'P^ 

ypi+---+p n =d p\ +—+p n =d-l 

+ •••+ 2 ci(pi,---,p B )r(0 (pi -"-'* ,) , </>l, i = 1,2,- ■• ,n.l. 
Pl+-+P„=l J 



We endow this space with the scalar product < •, • > 2 , given by 



n n 

= m p+q + ^ ^ c ei+q crijtl[ q=ej ), 
i=i ;=1 

W > 1, l9l > 1. 

Because 2 (A*i (*))*" • • • (AX B (s)) p " = 2 (AXiO)r ■ ■ ■ (AX B (j))M n ? jax,( s )*0), one clearly sees that 

0<j<< 0<s<t 

x Pl l x^ 2 1 ■ ■ ■ x%" 1 «-» yO>i.-.p») with |/?| > 1 and is an isometry between S \ and 52. An orthogonalization of 
{x^'xj 1 • • -x^pX^'xj 1 ■■■x n 1 ,-- - , x~ l , l,x\, ■ ■ ■ ,x n ,x\,x\X2, ■•■ ,x\, • ■ ■ }in Si gives an orthogonalization of 

ry(l,0,- ,0) ... y(0,-,0,l) y(2,0,",0) y(l,l,0,-,0) ... y(0,-,0,2) ...j 

It is well known that orthogonal polynomials of several variables are not unique (cf. Dunkl and Xu (2001)). In 
the remaining of the paper, {H p ,p e N"} is a set of pairwise strongly orthogonal martingales given by the previous 
orthogonalization of \Y P , p). It is also worth to emphasis that all deduction procedures and results are the same once 
the orthogonal martingales are determinatively given. 

3. Representation properties 

3.1. Representation of a power of a Levy process 
For notation simplicity, here and hereafter we set 

XP'V) = J] (*Xi(s)) Pt , Pi>2, i = 1 , 2, • • • , n. 

0<s<t 

and for convenience we put xf\t) = X,(f). Pi,qi(i = 1, 2, • • • ,ri) are all some nonnegative integers. Note that not 
necessarily X,(f) = 2 AX,(s) holds; it is only true in the bounded variation case with I = O. If S = O, clearly 

0<s<l 

[X h Xi](t) = Xf\t). The processes X^ = {Xf°'\t), t > 0}, p, ■ = 1, 2, • • • , are again Levy processes. They jump at the 
same points as the original Levy processes. 

We have E[Xi(t)] = E[Xf\t)] = fm i;1 < oo and by Protter (1990,p29), that 



E[Xf'\t)} = E 



^ (AXi) Pi = t I x Pi v(dx) = nii-p.t < oo, Pi > 2. 

Q<s<t J J 



Therefore, we can denote by 



yf'\t) d =l f X™(0 - E[xf\t)] = Xf\t) - m,, Pi t, Pi = 1 , 2, 3, • • • 

the compensated power jump process of order p„ and if' is also a normal martingale. 

We will express {X x {t + f ) - Xi(r )) 4 ' (X 2 (t + f ) - X 2 {to)f 2 ■ ■ -(X n (t + t ) - X n (t )) k ", t ,t > 0, k t = 1, 2, • ■ • , 
i — 1 , 2, • • • , n, as a sum of stochastic integrals with respect to the special processes Yf°'\f), i — 1 , • • • ,«,/?,■ = 1 , • • • , 

Using Ito's formula(Protter, 1990, p. 74, Theorem)we can write for fe,- > 2, i = 1, 2, • ■ • ,n, 

(X,(t + t ) - Xdto)) 1 " (X 2 (t + to) - X 2 (t Q )) k i ■ ■ ■ (X n (t + t ) - X n (to)) k » 

= 2 JC ki(Xi(s + to) - Xi{to)) k >- x U(Xj(s + t ) - Xj(to))dX,(s) 
i=i j& 

+3 2 C vlMki - l)(Xi(s + to) - Xi(to)) k -- 2 U(Xj(s + to) - Xj{to)) k 'ds 

i=l j*i 

+ 2 C "fMW + fo) - x i(to)) k '- l (Xj(s + to) - Xj(to)) k r l n (Xe(s + to) - Xdt )) k 'ds 

l<i<j<n 

+ £ {£ l(Xi(s + to) - Xi(t Q )) k < - (Xi((s + to)-) - Xi(t Q )) k q U(Xj(s + to) - Xj(t )) k > 

0<s<t (i=l j±i 

- 2 ki(Xi((s + to)-) - Xi{to)) k >- 1 U(Xj(s + to) - Xj(to)) k <AXi(s + t ) 
i=l j*i 



Z kiiXiiu) - Xjito)?'- 1 U(Xj(u) - Xj{t a ))dXf \u) 

i=l J '° j*i 



+iZ^-i) 



j£ «/ + fo) - Xi(<b))*- 2 + to) ~ Xj(t )) 



+ Z crfMj 

\<i<j<n 



t(Xi(t + t ) - Xi{to)) k '-\Xj{t + to) - Xj(to)) k r l n (Xe(t + to) - X ( (to)) kt 



- j£ sdUXiis + to) - XiiHdf'-HXjis + f ) - Xjtorfr-i ft (X ( (s + t Q ) - X { (t Q )) k 'ds 



+ Z { Z + fo)-) + AXKs + fo) - ^(fo)) t; - (X,-((s + fo)-) - X,-(<b))*»] 

xU(Xj(s + to)-Xj(to)) k ' 

- 2 + fo)-) - X^o))*' -1 U(Xj(s + to) - Xji^f 'AXiis + to) 



Z X o ° ktiXM - Xi(to)) ki - x U(Xj(u) - Xj{to))dXf\u) 



+2.2^(^-1) 



f(X,(f + fo) - X,(fo)) fc - 2 n(^(f + fo) - Xj(t )) k 



j*> 



- j£ «/ + fo) - Xiito)?'- 2 Y\(Xj(s + to) - Xj(t )) k 



+ _ z 4^ 



l</< j<n 



t(Xi(t + to) - Xiito^-HXjit + to) - Xj(to)) k r l n (Xe(t + fo) - X { (to)) ke 
- j£ sdUXiis + to) - X^f'-'iXjis + to) - Xj(t Q )) k '- 1 n (X((s + to) - X ( (t Q )) k 'ds 
+ Z Z Z ( 5 ) + f °)-) - n(^/(* + fo) - *,«*>))*' (AXj(j + fo)) f 



= Z C° +< - W)* i_1 U(Xj(u) - Xj(t Q ))dXl l \u) 



f(x,(f + f ) - x,a )f'- 2 nw< + fo) - Xj(t )) k 



J f ai + fo) - Xifo))*- 2 P(! + fo) - X/ft>)> 



+ Z ofMj 
\<i<j<n 



t{Xi(t + to) - X^f'-'iXjit + to) - Xjdo)?'- 1 n W + fo) - ^(fo))* f 



- C sdliXiis + to) - Xiito^'-HXjis + to) - Xj(to)) k r l U W + fo) - ^(fo))^) 



+ z zz ; (^wmw-wrnw-w 

t <u<t+t i=l 1=2 \ 1 I j*i 
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i=l(=2\ 1 J j*i 



i=l 



- JT' «/ + t ) - Xtto)?'- 2 U(Xj(s + t ) - Xj(t )) k 



+ 2 cr-Mj 

\<i<j<n 



t(Xi(t + to) - Xi(to)r- l (Xj(t + to) - Xj(to)f r l n (X((t + to) - X e {to)f 



J Q ' sdliXiis + to) - X^f'-'iXjis + to) - Xj(to)) k r l n (X e (s + to) - X ( (t Q )) k 'ds) 



(2) 



Lemma 1. The power of an increment of a Levy process, (X\(t + to) - Xi(to)) k ' (X 2 (t + to) - X 2 (to)) k2 ■ • ■ (X„(t + to) 
X„(to)) k ", has a representation of the form 

{X x {t + to) - X x {to)) h {Xi{t + to) - X 2 (t Q )t- ■ • • (X„(t + to) - X n (to)) k " 
= f {k Kt,to) 



dY, 



(^)-"rf^ , (^JT"- ,, (»4..i) 



■■■^ , * l) tt- 1 * 1 )-"rflf 1 ' ) ft- 1 ^lf'- ,) fe,l) 



,r/ ' f '' Y '' /,f 'C',-,, . :-:,..,.-.,„... 



(3) 

r2/in>*l + ''' + *»>\ 



(/, A)J ' " ' ' f ii,?,,;'" ; f i m; i^'" , ti„,q im ) are deterministic functions in L Z (W.+ m ), 
and k — {k\, ■ ■ ■ , k n ). In addition, the index m controls the number ofYj(t), i — 1, • • • , n, chosen from Y\(t), • • • , Y n (t). 
After the m is fixed, (z'i, • • • , i m ) indicates an arbitrary subset of the integer set ( 1 , 2, • • • , n). After the (z'i, • • ■ , i m ) 
is fixed, for ij e (z'i, • • • , i m ), the (f; i, • • • , fy, 9i .) wz'Z/z double index {ij, ■) indicates the time-points chosen from the 
time-points (t\, ■ ■ ■ , t^.) corresponding to (y,.(fi), • • • , F;.^.)). The meaning of power index Pi jt . is similar. 

Proof Representation (3) follows from (2), where we bring in the right compensations, i.e. we can write 



iil k ; \r\x i (s-)-xAto)) k -- e n(x j (s-)-x j (to)) k 'dxf\s) 

i=U=l\ <■ I '° j*i 

= 2 2 [ k l ) f + '° - Xi{to)) k -- e ll(Xj(s-) - Xj(to)) k >dYf(s) 

i=\(=l\ <■ l J '° i*i 



+ ZZ ;' \mul {X i (s-)-X i (to)) ki - t mx j (s-)-Xj{to)) ki ds 



= 2 2 

1=1 1=\ 



" b t ki \ rt+to 



n ki-1 / Z-. \ 

,?„?,(;)' 



I \Xiis-) - Xi(t )) k >- { U(Xj(s-) - X^f'dYfis) 



mtetiXiit + to) - Xi(t )r- e U(Xj(t + to) - Xj(t )pds 

-zz( 5 ) mu iT^^(^-)-^)) w nw-)-^/<b)) tj ) 

!=1 



(4) 
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Combining (2) and (4) gives 

n(x«(* +*>)-*,■(*>))*" 

= \ 2 ojM k i - !) W + fo) - X,(fo))^ 2 I! W + fo) - *,■(<(>)> 
i=i L m 

- JT' srf ((Xi(s + fo) - X«(fo))^ 2 n.W + to) - Xj(t )) 



+ 2 

l</< 7<« 



J*' 

/\,V> i- /(il - n. ; i/ (1 i/ 



?(*,.(? + f ) - x^f-^Xjit + to) - Xj(t )p- 1 n »(f + fo) - ^(f )r 



- j£ «/(*,■(* + fo) - Xi^))*'- 1 ^ + f ) - Xfayfr-i J] + fo) - X t {t )) h ds 



+ 2 i( *} )f\X i (s-)-X i (t Q )) k <- e U(Xj(s-)-X j (to)) k 'dYf ) ( S ) 

i=\(=l\ 1 I j*i 
n kj-1 I h. \ 



+ 22 ; maa*H)-%)rnw+ft)-w^ 

- 2 X? ( J ) m« «/ - Xfttifrt U(Xj(s-) - Xjito))^ 

n 

/=i 

(5) 

The last equation is in terms of powers of increments of X t which are strictly lower than kj. So by induction represen- 
tation (3) can be proved. □ 

Notice that taking the expectation in (3) yields 



y\(X i (t + t )-X i (t )) k < 



/ w (/,f )=/ (t) (f), f,f >0, 



which is independent of to- 

Moreover, it can easily be seen that 

f%l,L-,P lmi ;~; Plm ,L~,P lm , q , m ) ( - t ' fo; ^' ' ' ' ' f '-i*. ; ' ' ' ; fi - 1 ' ' ' ' ' Xi ^iO 
are just real multivariate polynomials of degree less than fc, and that we have 

€L-,P, imi ;-;p im . 1 ,-,A,„„ im ) (f ' fo; fi " 1 ' ' • ' ' ; ' ' ' ; fi »- 1 ' ' ' ' ' W = °' 

whenever + • • • + > fy. 

Because we can switch by a linear transformation from the F- Pi "" 2) (f;,„,2) to Hj Plm - (t imt2 ), it is clear that we also 
proved the next representation. 

Lemma 2. 77ze power of an increment of a Levy process, (X\(t + to) - X\{to)) ki {X 2 {t + to) - X 2 (to)) kl • • • (X„(t + to) - 
X„(to)) k ", has a representation of the form 

(Xi(t + to) - X 1 (to)) k '(X 2 (t + to) - X 2 {to)) h ■ ■ ■ (X„(t + t Q ) - X„(to)) k « 

\k\ 

= / w (f,fo)+2 2 C'° hf}{t,to\h)dm{h) 

\k\ 

+ 2 2 C £ h«l 2) (t, to; h , h )dm<t 2 )dw ( h) 

a=l p 1 +p 2 e!>!j 
+ ■■■ 



+ 2 2 I X 1 • • • X; H *£.-•*><'• ^o; fi, f 2 , • • • , f,*,)^(f W ) • • • dHnvxuvw 

a=\ pi+--+p\ k \eN d 
\k\ \k\ 

2 2 2 T" 1 f ' ' ' r™"' *£ ... , )< f > ^ f i> f 2' ■ ■ ■ • t m )dHP^t m ) ■ ■■dHPHt 2 )dHP'(t 1 ) 



=1 Pl+-+p m 
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where the index m controls the number of integral variables in each multiple integral, p t = (/?;,,••• ,PiJ € Ng, 
NJJ = {p e Ng : \p\ = d) and /i* j(f, fo; fi, ?2, • • • , f m ) aw deterministic functions in L 2 (M'"). <9f/zer notations are 
the same as Lemma 1. 

3.2. Representation of a square integrable random variable 

We first recall that {//. , i = 1, 2, • • • , n; p t = 1, 2, • • • } is a set of pairwise strongly orthogonal martingales, 
obtained by the orthogonalization procedure described at the end of Section 2. 

We denote by 

jf (!>,,■■■, p J = J/T g L 2 (Q) ; 

xr r • ■ • c f ° ; fi ' f2 ' ■ • • ' t ^ dHPm ^ ■ ■ ■ dH>*{t 2 )dHHh) (6) 

p°eN", /=l,2,---"m) 
We say that two multi-indexes 

(Pu---,Pm) and (pu---,Pm) 
are different if m + m or when m = m, if there exists a subindex 1 <l<m = m, such that /?/ ^ and denote this by 

(Pi,- - ,Pm) * (P\,--- ,Pm) 

Proposition 1. If 

(Pi,-- - ,Pm) * (PW- ,Pm) 

then 

jff>(Pl,~ ,Pm) j_jtff(Pl,~ ,Pm) 

Proof Suppose we have two random variables K e Jif ipi '' "> Pm> and L e j/jrffo"- <P*> . We need to prove that if 

(Pu- ■ ■ ,Pm) + (pu ■ ■ • ,Pm) 

then K±L. 

For the case m = in, we use induction on m. Take first m-m-\ and assume the following representations for K 
and L: 

poa poo 

K= f(t x )dH Pi ( tl ), L= g{t x )dH p \h) 
Jo Jo 

where we must have p\ + p\, By construction H Pl and H pl are strongly orthogonal martingales. Using the fact that 
stochastic integrals with respect to strongly orthogonal martingales are again strongly orthogonal (Protter, 1990,Lemma 
2 and and Theorem 35, p. 149) and thus also weakly orthogonal, it immediately follows that K±L. 

Suppose the theorem holds for all 1 < m = m < n — 1 . We are going to prove the theorem for m = m — n. Assume 
the following representations: 

K = C'° £ •'• fj-*'*'* 1 '* 2 '"' ,tm)dH p *{t m )---dH p i{t 2 )dm{h) 

= ^a{h)dm{tl) 

L = £'° £ ■ ■ ■ g ( £,.., p Jt, t ; tuh,— , t m )dHHt m ) ■ ■■dHHti)dH p '(t l ) 
= /3(h)dHHh) 

There are two possibilities: (1) pi = p\ and (2) p\ + p\. In the former case we must have that 

(P2, ■ ■ ■ ,Pm) * (Pi, ■ ■ ■ ,Pm) 
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and thus by induction a(t\)±B(t\), so that 

E[KL] = E\ \ a s B s d <H P, ,H P ' > 



f 

Jo 



E(a s B s )d < H Pl ,H Pl > s =0. 



o 

In the latter case we use again the fact that stochastic integrals with respect to strongly orthogonal martingales are 
again strongly orthogonal (Protter,1990,Lemma 2 and and Theorem 35, p. 149) and thus also weakly orthogonal. So 
it immediately follows that K±L. 

For the case m + m, a similar argument can be used together with the fact that all elements of every J^to.-.^ 
£ > 1, have mean zero and thus are orthogonal w.r.t. the constants. □ 

Proposition 2. Let 

In m | 
Y\ Y\(Xi(tj) - X,(ii ,)) 4 : m > 0,0 < t < h < h < ■ ■ ■ < t m , * u , ■ ■ ,£„, m > 1 J 

then we have that & is a total family in L 2 (Q, i.e. the linear subspace spanned by & is dense in L 2 (Q, 

Proof Let Z e L 2 (Q., J£") and Z±ZP. For any given s > 0, there exists a finite set {0 < s\ < ■ ■ ■ < s m ] and a square 
integrable random variable Z e € L 2 (Q, cr(X(si), X(s2), • • • , X(s m ))) such that 

E[\\Z-Z s \\ 2 ]<s. 

So there exists a Borel function / such that 

Z s = f e (X(si), X(s 2 ) - X( Sl ), • • • , X(s m ) - X{s m ^)). 

Because the polynomials are dense in L 2 (M",P o X(t)~ l ) for each t > 0, we can approximate Z B by polynomials. 
Furthermore because Z±0*, we have E[ZZ e ] = 0. Then 

E [\\Z\\ 2 ] = E[Z ■ (Z - Z e )] < ^E[\\Z\\ 2 ]E[\\Z-Z S \\ 2 ] < ^sE[\\Z\\ 2 ], 

and Letting s — > yields Z = a.s. Thus & is a total family in L 2 (Q, □ 
We are now in a position to prove our main theorem. 

Theorem 1. ( Chaotic representation property ( CRP)). Every random variable F in L 2 (£2, has a representation of 
the form 

f= e(f )+ z S 2 rf - 

m=l J=l p]H — Hp„,eNJ 

Jo"' fSl-fJt' fo; f1 ' f2 ' • • • ' tm)dH p "(t m ) ■ ■ ■ dm{t 2 )dm{h) 

where the /"*... kh h,h,--- , t m )'s are functions in L 2 (M+). 

Proof Because & is a total family in L 2 (f2, ^"), it is sufficient to prove that every element of & has a representation 

n m 

of the desired form. This follows from the fact that & is build up from terms of the from \\ \\(Xi(tj) - 

i=\ y=l 

wherein every term has on its turn a representation of the form (6), and we can nicely combine two terms in the 
desired representation. Indeed, we have for all ki, /,■ > 1, i = 1, • • • ,n, and 0<f<s<M<v, that the product of 

n 

Y[(Xi(s) - Xi(t)) ki (Xi(v) - Xiiu)) 1 ' is a sum of products of the form AB where 

i=\ 

A = f I,'" ■ ■ ■ t~ l ~ h ™,-,p m P> *> ^ t2 > ■ ■ ■ • ^dHO-itm) ■ ■■dHK{t 1 )dm{h) 
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and 



B = rr.-p-'-/,' 1 ' 

Ju Ju Ju 

where m and m are two integers. 
We can write 



A) (v, m; mi, M2, ■ ■ ■ , Ur h )dHP"{t ih ) ■ ■ ■ dHP\t 2 )dm(h) 



ab = £ r~ • • • I"* -1- r r ~ • • • x' m ~ i_ /,; » ... «; ^ «2, ■ • ■ , ««) 



. . . i ; j 

f; fl ' f2 ' ' • ' ' t m) dHPm <Jm) ■ •"• dH^{t 2 )dH^{h)dHHtih) ■ ■ ■ dHHh)dHHh) 

= XX' '"X"" X'" X' '"X"' 1 hu,v]( u l)hu,m](U2) ■ ■ ■ lo^lOm) 

l(M](*l)l(«i](*2) • • • l(W„-i](fm)^... p nl) (v, M; U\,U 2 , ' • • , u m)h^ w .. pmj {s, f, t\, t 2 , ■ ■ ■ , t m ) 

dm- (t m ) ■ ■ ■ dH P2 {ti)dm (tx)dHHtin) ■ ■ ■ dm (t 2 )dH pl (h) 
and the desired representation follows. □ 

Theorem 2. (Predictable representation property (PRP)). Every random variable F in L?{Q., has a representation 
of the form 

OO 

F = E(F) + Z Z X OP(s)dHP(t m )(s) 

where <$ p (s) is predictable. 

Proof From the above theorem, we know that F has a representation of the form 
F - E(F) 

oo oo _ . _ 

= zz z X X ■• ■j[^ 1 "y£t..^i.»2,---,»m)rfff fc (»m)---rfff*fe)rfff' i (»i) 

m= 1 fif=l /»i+---+/> m €N2 m 



2 Z X 

i=0 P2+-+/» m €N! 



z z £°f ( £ ) (h)dH^(t l )+^ z X 

Jo™" 1 , )('i-'2. ■ ■ ■ • t m )dH?>"(t m ) ■ ■ -dH^h)] dH*(ti) 

I, ..I, JT i«<") 



+ Z Z X X /( h-» ,tm)dHP"{t m )---dH»{t 2 ) 

k=0 p2+-+Pm^l ' 

oo 

z z X visW'^xs) 



dHP'(h) 



which is exactly of the form we want. □ 

Remark 2. Because we can identify every martingale M e fy 1 with its terminal value M M € L 2 (Q., and because 
M, = E[Moo|^i], we have the predictable representation 

which is a sum of strongly orthogonal martingales. 

Another consequence of the chaotic representation property, is the following theorem: 
Theorem 3. We have the following space decomposition: 



L 2 <A JT) = 



\d=\ peW. 
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Remark 3. The Levy-Khintchine formula has a simpler expression when the sample paths of the related Levy process 
have bounded variation on every compact time interval a.s. It is well known(Bertoin, 1996,p.l5), that a Levy process 
has bounded variation if and only if 2 = 0, and J(l A ||*||)v((f;t) < oo. In that case the characteristic exponent can be 
re-expressed as 

4/{0) = id-0 + \ (exp(i0 ■ x) - 1) v(dx). 

Furthermore, we can write 

Xj(t) =dt+ Yj aX,-(j), t > 0, i = 1, 2, • • • , n. (7) 

0<s<t 

and the calculations simplify somewhat because S = and for k > 1, 

2 Jf fcQftj + f ) - XKfo)) 4 '- 1 Il^i + fo) - X/f )¥W 

i=l 

= Z ( Z + <*>)-) - ^(fo))^ 1 1! W + fo) - Xj(t )) k ' AXiis + t )\ 

0<s<t (i=l j±i ) 



4. Examples 

Multidimensional models with jumps are more difficult to construct than one-dimensional ones. A simple method 
to introduce jumps into a multidimensional model is to take a multivariate Brownian motion and time change it 
with a univariate subordinator (refer to Cont and Tankov (2004)). The multidimensional versions of the models 
include variance gamma, normal inverse Gaussian and generalized hyperbolic processes. The principal advantage 
of this method is its simplicity and analytic tractability; in particular, processes of this type are easy to simulate. 
Another method to introduce jumps into a multidimensional model is so-called method of Levy copulas proposed by 
Kallsen and Tankov (2006). The principle advantage in this way lies in that the dependence among components of 
the multidimensional Levy processes can be completely characterized with a Levy copula. This allows us to give a 
systematic method to construct multidimensional Levy processes with specified dependence. 

In the following first and third examples, we define a multivariate gamma process and a multivariate Meixner 
process by using Levy copulas, and furthermore discuss their orthogonalization procedures. All the concepts and 
notations are adopted from the Kallsen and Tankov (2006). In particular, for n > 2, the Levy copula F{u\, ■ ■ ■ , u„) : 
E" -» R is taken as 



F{u\,--- ,u„) = 2 2 



(^{«i -«„>0} - (1 - 7)'{ki"«„<0))- (8) 



It defines a two parameter family of Levy copulas which resembles the Clayton family of ordinary copulas. It is in 
fact a Levy copula homogeneous of order 1, for any 9 > and any rj e [0,1]. 

In addition, we know that if the tail integrals f/,(x,), i = 1, • • • ,n, are absolutely continuous, we can compute the 
Levy density of the Levy copula process by differentiation as follows: 

v(dx u --- ,dx n ) = d\ ■ ■ ■ d n F\ il=Ul(xlX ...£ n=UnUn) vi(dx\) ■ ■ - v\(dx n ) (9) 

where v\(dxi), • • • , v n (x n ) are marginal Levy densities. 



4.1. The multivariate gamma process 

In the literature, the multivariate gamma distributions on W have several non-equivalent definitions(refer to John- 
son and Balakrishnan (1997)). Here we consider only a multivariate gamma process by using copula. The multivariate 
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Gamma process G(t) = (Gi(f), ■ • ■ > G„(f)) r is a multivariate Levy process with the marginal distribution density 
functions of G,(/)> 1 ■ = 1.2, ■• ■ , n given by 

fG,(t)(Xi) = —^—Af'xf'^expi-AiXi}, 

r(r,o 

x,- > 0, /I,-, y,- > i = 1, 2, • • ■ , n. 
The corresponding marginal characteristic functions are given by 



i = 1,2, • • • , n. 



The corresponding marginal Levy measures are given by 

Vi{dXi) = ^exp{-AiXi}I ( Q t oc,){xi), i = 1, 2, • • • ,n 
The corresponding Levy measure is given by 

Fir; « 

v/»(Jx) = <9i ■ • ■ 5 B F| fl=£/l(xi)j ... A=£/i>ta) ^-ex/>{- 2 /tiXi}/ ( o,oo)»(*i, ■ • • , 

1=1 

i=l 

where <ijc = dx\ ■ ■ ■ dx„. The n-dimensional Gamma processes are used i.a. in insurance mathematics(Dickson and 
Waters, 1993, 1996; Dufresne and Gerber, 1993;Dufresne et al., 1991). 
We denote by 



Gf\t) = £ (AGi(i)f, Pi > 1, i = 1,2, • ■ • ,n 



0<s<( 



the power jump processes of G,(0- In addition, set G (pw "> p "\t) = (G\ Pi) , ■ ■ ■ ,G^ n) ), where (p u - ■■ ,p„) e N". Using 
the exponential formula (Bertoin, 1 996), and the change of the variable z = (x^ 1 , • • • , x p "), we obtain for p 1 H h/?„ > 1 



E[exp{ie T G (p ^' p "\t))] 
= exp 



= exp 



I n \ nn n 

t L \exp(i 2 OiX P ') - 1 di • ■ ■ d n F\ il=UlUll ...^ n=UM '^—exp{- £ AiX^dx 

v + V i=\ I n*i 1=1 J 

n« » ] 

tj n {exp(\e T z) - 1)5! • ..5 n F| fi=£/ife i /P1) ... fn=£/dfe i /P „ ) -^ e xp{- g Ajz/^Jdz 
which means that the Levy measure of G {pu '" ' Pn) is 

[Ir- 

i=l 

Introduce power jump processes of the form 

def 

0<.5<f 

and then define the Teugels martingale monomial 

G(t) (pi '" ' p " ) = f G(f) (/ "'' , ''"" ) - E[G(t) (p "-' p " } ] = G(t) (p, '-' p " } - m p t. 
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G(f) (P„ ■,/>„) <¥ 2 (AGi(s))"> • • • (AG n (s)) p 



Because 



2 (AGi(s)) Pl ■ ■ ■ (AG n (s)) p " 

Q<s<t 



= f Jr- x i l ■ ■ ■ x " dl ' ' ' 5 »- F lft=f/ 1 (x 1 ),-,f„=c/„(x„) i T-ex;'{- 2 i,x,}c/x, |/»| > 1, 
Next, we orthogonalize the set G p of martingales. So we are looking for a set of martingales 



H" 



(10) 



such that // p is strongly orthogonal to H p , for p + p. 

The first space S \ in the gamma case is defined as follows 

(d 

S\ = \Yj^j c k{pu- - ,p n )x Pl ■■■x p n " +cq+ ^ ■ - ■ , i n )x'l 

(*=1 \p\=k (ii,-,/„)e{0,-l}»,|i|>-(B-l) 

d e {1,2, • • ■ },c,Oi, •••,/?„) e K,j = -1,0, ••• > 0,/ = 1, ••• ,n;i = 0'i, •• • ,/„)} 
which is endowed with a scalar product < •, • >i, given by 
<P(x),2(x)> 1 

p+oo r+co 



= X °° " ' i) °° p ( x )2( x )<9i • • • 5„F| fl=£/l fe),-,f„=£/„fe) U(nxdexp{- 2 i;Xi)c/x. 



i=l 



Note that 



1+91-1 v P„+q„-l; 



= Lj^i 1 ^ 1 ■■■ x n' +q " d i ■■■dnF\fr= Ul (xi),---&=u n (x n )\nyijexp[- t^iXi}dx 
\p\,\q\>l- 

Thus we can construct the other space S 2 which is the space of all linear transformations of the Teugels martingale 
monomials of the multi-dimensional Gamma process, i.e. 



S 2 = 



J] a d {pu--- ,p n )G(t) (pi '-' p » ) + J] aa-iQn,--- ,p n 

\pi +—+p n =d pi+-~+Pr,=d-\ 

+ •••+ J] alto,--- , Pn )G(t)^"- p »\ d>i\. 

Pl+-+P„=\ j 



)G(t) {Pi '- 



(j>l,~,Pn) 



= E 
= E 



endowed with the scalar product < •, • >2, given by 

< Q(pi,",P») t Q(qi,",q«) > 2 

[G(Pi,-,P») 5 G(9i,-,9»)j(l)j 

r+1+1 I n \ n 

= Jr» x f +? '~ • • • 3i ■ • ■ 3 n F| fl=£/lfa) ,..., fn=£/nfe) I n 7,1 2 ^a ( }Jac 

So one clearly sees that x p '~ l x P2 ~ l ■ ■ ■ x p "~ l «-» G^ Pl '"' ,p "^ is an isometry between 5 1 and An orthogonalization of 
{x~[ l x 2 l ■ ■ ■x~l l ,x^ 1 x^ i ■•■x n 1 ,-- - , x~ l , l,x\, ■ ■ ■ ,x n ,x\,x\X2, ■■■ ,x\, • ■ ■} in S\ can give the multivariate polynomi- 
als, so by isometry we alsocan find an orthogonalization of {G< 10 ■ °), . . . j (j(0,-,o,i) ) G (2,o, -,0^ 3(1,1,0, -,0^ . . . . . . } 
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4.2. The negative multinomial processes 

The next process of bounded variation we look at is the negative multinomial processes, sometimes also called 
Pascal processes. Here the conception of negative multinomial processes can be found in Johnson et al.(1997), and P. 
Bernardoff(2003). 

We define a negative multinomial distribution on N!!. Its distribution is £ P r °bk5k, where probkSk denotes the 
probability measure concentrated at k = \k\,ki, • ■ ■ , k„], 

r(r + £ kd 

prob k = P(fci, ••-,£„) = - - , l= ] — A'UQiAif, 

kx\k 2 \ ■ ■ ■ k n \Y{f) ij 

ki = 0,1,2,-- • , ? = 1,2, - ■ . 

where < A < 1, < pAj < 1 for i = 1,2, ■ • • ,n and A + p{A\ H 1- A n ) = 1. 

This type of n-dimensional levy processes P = {P(t), t > 0} where P(t) = (Pi(t), Pi(t), ••• , P n (t)), has a character- 
istic function given by 

E [expiiff ■ P(t))] = Yj Probke^ 01 ■ ■ ■ e ikA 
keN'< 



\1 -pAAie^ + ■ ■ ■ + A„e ie ») ) 
The corresponding Levy measure v(k\ , • • • , k n ) is given by 

v(h, ■ ■■,k n )= ^ - ] } \ fj (M)*' , 1*1 = f h + ■ ■ ■ + kn. 

Let us denote with 

G(t) (pi '" ' Pn) d = J] (a^K*))" 1 " --(APnWr, |J>| > 1. 

0<.s<r 

the power jump processes of P and with gffi.-.^) = .-•/>«)(f) ? f > 0) the corresponding processes of Teugels 
martingale monomials. 

We look for the orthogonalization of the set [Q(pi"--P"\ \p\ > 1} of martingales. The space S\ is now defined as 
follows 

Si = \J]Yj c k(Pu-- ,p n )k Pl ■■■k Pn +c + 2 ,i n )k l l ■■■k i ;\ 

\k=\ \p\=k (i 1 ,-,«€{0,-l)»,|/|>-(n-l) 

d e {1,2, ■ • -},c/(/?i, • • • ,p„) e R,j = -1,0, •• • , d; k t e N, i = 1, • • • ,n;i = (/[,-• • ,/„)} 
endowed with a scalar product < ■, ■ >i, given by 

< p(k),R(k) > 1= V p(k)R{k)h • ■ -^rr-rr fj (MO* ■ 

Note that 

<*f- 1 ...c _1 ,*f- 1 -c- 1 >i = zc^-c^'^nw) 1 ' 

t£NJ "' 1=1 

> 1, ;', j = 1,2, ••• ,n 
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The other space S 2 is the space of all linear transformations of the Teugels martingales of the negative multinomial 
processes, i.e. 

Si = J J] a d(PU--- ,Pn)Q(t) iPl '"' P " ) + Y a d-l(Pl'"-,Pn)Q(f) iPl -"'- P " ) 
I p\-\ — \-p n =d pi+--+p n =d-\ 

+ •••+ Y adpu--- ,Pn)Q(t) (pi '-' p "\ 
pi+-+p„=i 

deN, (pi,---,p„)e N", aj (p u ■■•,/>„) e M, 7 = 1,2, • • • , d}. 
and is endowed with the scalar product < •, • >2, given by 



keN- 



So one clearly sees that k P l 1 ' l k P2 ' 1 ■ ■ ■ k^ 1 «-> Q(Pi" m -P*> is an isometry between 5 1 and 5 2- An orthogonalization of 
{k\ l k 2 l ■ ■ ■ k~\ v k[ l k^ ■ ■ ■ k n • • • .fc" 1 , l.fci, • • • ,k n ,k 2 , k\ki, ■ • ■ , k\, ■ ■ ■ } in S2 gives the multivariate Meixner poly- 
nomials(Griffiths(1975), Griffiths and Span6(2008), and Koekoek and Swarttouw(1998)), so by isometry we also find 
an orthogonalization of 

{g(l,0,-,0) ) . . . ) q(0,-,0,1) j q(2,0,-,0) j q(1,1,0,-,0) ) . . . > G (0,- -,0,2)^ . _ . j 



4.3. r/ze Multivariate Meixner process 

A multivariate Meixner process M(t) = {M\{t), Mi(i), ■ ■ ■ ,M n (t)) T ,t > is a bounded variation Levy process 
based on the infinitely divisible distribution. We use the copula to construct a multivariate Meixner process. Here the 
marginal density functions are given by 

f I \ (2cos(a i /2)) 2m i , s.r-y , • n,2 

m, ad = 2 ^r(2m') expiaiXiWinii + , 
Xj e (-00, +00), i= 1,2, ■ ■ ■ ,n 

The corresponding distribution is the measure of orthogonality of the Meixner- Pollaczek polynomials (Koekoek and 
Swarttouw, 1998). The Meixner process was introduced in Schoutens and Teugels (1998). In Grigelions (1998), it 
is proposed for a model for risky assets and an analogue of the famous Black and Scholes formula in mathematical 
finance was established. The marginal characteristic functions of M,(f), ;' = 1, 2, • • • , n, are given by 

E[exp(6 i M i (t))] = { C h Z (ai/2 \nT ' ' i = ^'-> n - 
\cosh{{9i - ia,)/2) / 

In according to the results in Schoutens and Teugels (1998) and Schoutens(1999) and applying (11), its Levy measure 
can be calculated as: 



1 j ^ q n rii fl niiexpiaiXi) 



i=l 



= di---d n F\ h=Ul(xiX ... 4n=UM Y\mi\T{\ +ix i )\ 2 ^^^-dx i . 
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Also note that 



f~[ x 2 v(<ix) = 5i • • •3„F| fl=f / lta ),..., ft - f / n(Xn ) ]~~[ m,-|r(l + ix,)| 2 -^-dxi. 



Being completely similar as in the above two examples, we can orthogonalize the multivariate Teugels martingales 
for the multivariate Meixner process by isometry. 

Acknowledgment. The author thank a kind proposal given by Professor David Nualart for the initial version of this 
paper. 
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